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^ . Hadronic Z decay data taken with the ALEPH detector at LEPl are used to measure the three-jet 
Ph! rate as well as moments of various event-shape variables. The ratios of the observables obtained 
^ from 6-tagged events and from an inclusive sample are determined. The mass of the b quark is 
extracted from a fit to the measured ratios using a next-to-leading order prediction including mass 
effects. Taking the first moment of the distribution, which is the observable with the smallest 
hadronization corrections and systematic uncertainties, the result is 



mfe(Mz) = [3.27 ± 0.22(stat) ± 0.22(exp) ± 0.38(had) ± 0.16(theo)] GeV/c^ . 

The measured ratio is alternatively employed to test the flavour independence of the strong 
coupling constant for b and light quarks. 
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1 Introduction 



The advent of next-to-leading order (NLO) calculations for e"^e~ annihilation into quark pairs, 
which take full account of quark mass effects |Q, ||, ^, has opened up the possibility for further 
investigations of QCD, such as measurements of the flavour independence of the strong coupling 
constant ||, or measurements of the running 6-quark mass at the Z mass scale f^, ||. The 
6-quark mass is one of the fundamental parameters of the Standard Model Lagrangian, and can 
in principle be viewed as a parameter similar to the coupling constant. The renormalization 
procedure leads to the definition of a running quantity, the renormalized 6-quark mass, which is a 
function of the renormalization scale. An interpretation as a particle mass is difficult because of 
the fact that quarks are not asymptotically free states. Nevertheless, a definition of the mass as 
the pole of the quark propagator is used frequently, which, depending on the size of long-distance 
QCD effects, can be interpreted as the mass of an (almost) free particle. 

Most of the 6-quark mass determinations have been performed at rather low scales 0. It is 
therefore interesting to measure this parameter at a large scale. Previous measurements have used 
the ratio of three-jet rates in b- and wrfs-quark decays of the Z boson. For a running 6-quark mass 
of about 3 GeV/c^ a 3% deviation from unity is observed in this ratio This is because the large 
6-quark mass suppresses gluon radiation in a similar way to the suppression of bremsstrahlung for 
muons compared to electrons. 

In this analysis a set of event-shape observables has been employed in addition to the three-jet 
rate in order to study the 6-quark mass dependence of the measured ratios of the observables 
in b and uds initiated events. The suppression of gluon radiation by the 6-quark mass affects 
the event-shape distributions, and these variables have rather different behaviour with respect 
to the hadronization and next-to-leading order effects. The final result is derived using the first 
moment of the ys distribution, which is the least affected by hadronization and next-to-leading 
order corrections. 

In the following section the analysis method is outlined, then a description of the ALEPH 
detector is given, followed by a summary of the data analysis. In Section || the results for the 
measurements of the ratios of observables in b and uds events are given, which are then used to 
determine the 6-quark mass in Section |^. A test of the flavour independence of the strong coupling 
constant is described in Section |], and the conclusions are given in Section 0. 



2 The analysis method 

The running 6-quark mass enters into the theoretical prediction for the ratio 

^r = ^ , (1) 

where Ob and Oi are infrared and collinear safe observables at the perturbative level (quarks and 
gluons) for Z decays into b and light (/ = uds) quark pairs, respectively. The measured ratio R^^'^^ 
of the observables in 6-tagged events and in an all-flavour inclusive sample can be related to the 
quantities at the parton level via the following formula: 

i?^;^^ = {ObHbDbnVb + OcH,D,T,Vc + OiHiDiTiVi) / O^""^ . (2) 
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Here if^;, D^, are the corrections due to hadronization, detector effects and tagging, respectively, 
and Vx are the purities of the tagged sample, where x is the true flavour. The tagging corrections 
take into account biases introduced by the flavour tag. The denominator O^'^^^ can be rewritten 

as 

O^"^^'^^ = Rb Ob Ht Dt + Rc Oe D,+ {l~Rb- Rc) Oi Hi Di , (3) 

where Rb(c) is the ratio of the partial width of the Z to 6(c) quarks and the total hadronic width. 

Taking all these ingredients, the relation between the measured ratio R^'^^^ and the ratio of interest 
rt . ^. ^ 

u 



Rli'^^ is given by 



Rb Rl^^^ Hb/i Db/i + Rc R^i"^^ Hc/i Dc/l + (1 — -Rft — -Rc) 



where Hb/i = Hi,/ Hi, D^/i = Db/Di, and the ratio -R^r* is set to 1. A small deviation of this ratio 
from unity because of the c-quark mass is considered in the systematic studies. All the correction 
factors and purities are obtained from Monte Carlo (MC) simulations. Because mainly ratios of 
corrections are involved, some systematic uncertainties cancel. 

Rli'^^ is extracted from the relationship @j and finally corrected for the contribution of 
anomalous triangle diagrams |T0[, in order to relate Rfl^^ to R^^^^ for which the perturbative 



calculations have been performed. The triangle contributions cancel out in the second ratio. 
However, they give a contribution of the order of 0.2% to the first. 

The following observables have been considered: 



The rate of three-jet events, where the jets are defined by the DURHAM clustering algorithm 
with the E recombination scheme [0. As shown in Ref. 0, the optimal choice of the 



resolution parameter |/cut for the mass determination is 0.02, because it represents a good 
compromise between the actual size of the mass effect to be measured and the size of 
backgrounds from other jet topologies. 

The first and second moments of the event-shape variables thrust T, the C parameter, the 
three-to-two jets transition value and the total and wide jet broadening, Bt and Bw 



The definitions of the observables can be found, e.g., in |]T2|, |T3[ and references therein. 
Moments have been chosen instead of distributions because of the better statistical accuracy 
with which the perturbative predictions can be evaluated by MC integration. Even higher 
moments would give smaller hadronization uncertainties. However, the NLO corrections 
typically get larger and the overall mass effects fall below the 3% level. 



A list of the leading (LO) and next-to-leading order (NLO) contributions to 1 — R^^"^ is 
given in Table |l| for all variables, both for the running and the pole mass schemes. They have 
been evaluated for a 6-quark mass of 3 GeV/c^ in the former and 5 GeV/c^ in the latter case. 
The computation of these terms is described in detail in Section |5.1| . It is found that for some 
observables such as the jet broadening variables the mass effect is rather large. However, the 
NLO corrections can also be sizeable, as in the case of the second moment of thrust and the first 
moment of the wide jet broadening. For such variables it would definitely be necessary to also 
compute the NNLO contributions in order to obtain a reliable perturbative prediction. Because of 
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Table 1: Leading (LO) and next-to-leading (NLO) order contributions to 1 — i?^™*, for the running mass (run) 
and the pole mass (pol) schemes. The contributions are evaluated for a running (pole) mass of 3 (5) GeV/c^. The 
strong coupling as{Mz) is set to 0.119. The values are given for the three-jet rate (DURHAM scheme) and for the 
event-shape variables thrust T, C parameter, the transition value 1/3 for three to two jets (DURHAM scheme), and 
the total and wide jet broadenings {Bt and Bw)- The indices indicate the first or second moment of the event 
shape variable. 






LO (run) 


NLO (run) 


LO (pol) 


NLO (pol) 


R3 


0.020 


0.010 


0.056 


-0.008 


Ti 


0.036 


0.019 


0.076 


-0.007 


T2 


0.017 


0.032 


0.043 


0.036 


Ci 


0.044 


0.022 


0.091 


-0.011 


C2 


0.021 


0.039 


0.052 


0.043 


l/3i 


0.032 


0.007 


0.071 


-0.021 


1/32 


0.015 


0.003 


0.032 


-0.007 


Bti 


0.117 


0.006 


0.188 


-0.074 


Bt2 


0.036 


0.112 


0.080 


0.123 


Bwi 


0.117 


-0.085 


0.188 


-0.183 


Bw2 


0.036 


0.016 


0.080 


-0.013 



these observations in the following only those variables are considered for which in both schemes 
the NLO contribution is clearly smaller than the LO term. This requirement selects the three-jet 
rate (-R3), the first moments of thrust (Ti), C parameter (Ci), the three-to-two jets transition 
value (1/3 J and the total jet broadening (-BtJ, as well as the second moment of the three-to-two 
jets transition value (2/32) and the wide jet broadening [By/^). 

3 The ALEPH detector 

The ALEPH detector is described in detail elsewhere |1l4| , [T5| . Here only a description of the 
tracking detectors is given, being the relevant ones for this analysis. Briefly, at the core of the 
tracking system is a silicon strip vertex detector (VDET). This has two layers, at average radii 
of 6.5 and 11.3 cm, each providing measurements in both the r-cf) and r-z views. The spatial 
resolution for r-cj) coordinates is 12 /xm for normal incidence and varies between 12 and 22 /im 
for z coordinates, depending on the track polar angle. The angular coverage of the VDET is 
|cos^^| < 0.85 for the inner layer and |cos0| < 0.69 for the outer layer. The VDET lies within a 
cylindrical drift chamber (ITC), which measures up to eight coordinates per track in the r-0 view, 
with a resolution of 150 /im. The ITC is in turn enclosed in a large time projection chamber (TPC), 
lying between radii of 30 and 180 cm. The TPC provides up to 21 three-dimensional coordinates 
per track, with resolutions in the r-0 and r-z views of 180 /xm and 500 /im, respectively. The 
three tracking detectors are surrounded by a superconducting solenoid producing a magnetic field 
of 1.5 T. 

For charged tracks with two VDET coordinates, a transverse momentum resolution of 
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^Pt/pt = 6 X 10 ^pr©0.005 {pt in GeV/c) is achieved. The three-dimensional impact parameter 
resolution is (25 + 95/p) /xm {p in GeV/c). 

Recently the LEPl data set has been reprocessed using improved reconstruction algorithms. 
In particular, a new VDET pattern recognition algorithm allows groups of several nearby tracks, 
which may share common hits, to be analyzed, in order to find the hit-to-track assignment which 
minimizes the track fitting for the event as a whole. The improvement on the hit association 
efficiency is more than 2%. Information from the TPC wires, in addition to that obtained from 
the pads, is used to improve the coordinate resolution by a factor of two in z, and by 30% in r-0 
for low momentum tracks. 

The tracking system is surrounded by electromagnetic and hadronic calorimeters, and together 
they are used to measure the neutral and charged energy ffow. 



4 Data analysis 
4.1 Event selection 

In this analysis data taken at the peak of the Z resonance from 1991 to 1995 are used. A standard 
hadronic event selection is applied, which is based on charged particles. A cut |cos^^t| < 0.7 



is imposed, where 9t is the polar angle of the thrust axis, computed from all charged and neutral 
particles as obtained from the energy-ffow algorithm |16[. This requirement ensures that the 
events are well contained within the VDET acceptance. According to the MC simulation the 
event selection is 61.7% efficient. Non-hadronic background, which is dominated by r+r~ events, 
amounts to 0.3% of this sample. After the selection, a sample of 2.3 million hadronic events 
remains for further analysis. 

The analysis also uses 8.8 million simulated hadronic events produced with a generator based 



on the JETSET 7 . 4 parton shower model fT^j . The production rates, decay modes and lifetimes of 
heavy hadrons are adjusted to agree with recent measurements, while heavy quarks are fragmented 
using the Peterson et al. model [|l8l. Detector effects are simulated using the GEANT package |jl9| . 



The MC events are reweighted in order to reproduce the measured values for the gluon splitting 
rates into cc and hh pairs, which are g^c = 0.0319 ± 0.0046 and gf,i = 0.00251 ± 0.00063 l2|. 



The observables described in Section ^ are computed using only charged tracks. This choice is 
preferred over taking both charged and neutral energy-ffow objects, as for this analysis it allows the 
reduction of the detector related systematic uncertainties to below the 1% level, without reducing 
the statistics. 

For the computation of the three-jet rate, an additional cut on the minimum energy of a jet of 
5 GeV is applied, which removes about 5% of all three-jet events. The three-jet rates are 17.7% 
(18.2%) in data (MC) for 6-tagged events, and 19.3% (20.4%) in the inclusive sample. The excess 
of three-jet events in the MC simulation had been observed previously |T2 |. 
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Figure 1: Distribution of the 6-tag variable for data and Monte Carlo b, c and uds events. The Monte Carlo has 
been normalized to the same number of events as the data. The last bin includes overflow entries. The insert at 
the bottom shows the ratio of data over Monte Carlo. 



4.2 The 6-tag algorithm 

The 6-tag algorithm, i.e., the criteria used to select b events, follows rather closely the approach 



described in detail in Ref. [^. Briefly, the presence of b hadrons is detected using a tag based 



on the long b hadron lifetime and the precision of the VDET. In contrast to [^], where the tag 
is applied separately to both hemispheres of the event, here the algorithm uses all tracks of the 
event. 

The actual selection of b events is obtained from a cut on the distribution of the confidence level 
Puds that all tracks of the event come from the main vertex. This distribution as measured in data 
and in MC is shown in Fig. |I], together with the expected contributions from different flavours. 
Very good agreement between data and MC is observed. The requirement — logj^g Puds > 2.2 
selects b events with an efficiency of = 80.5% and a purity of Vb = 83.1%. The backgrounds 
amount to Vc = 13.5% and Vi = 3.4%. 

With this selection, about 488 000 6-tagged events are found in the data before the three-jet 
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requirement. Out of about 435 000 three-jet events, 86 000 are tagged as bb. 



4.3 Detector and tagging corrections 

The observables must be corrected for detector effects, such as acceptance and resolution. This 
is achieved by computing the observables before and after detector simulation, imposing the 
same track and event selection cuts as for the data, except for the flavour tag. The ratios of 
the observables are computed for each true flavour, i.e., the quark flavour of the primary decay 
products of the Z boson. When computing the observables at the true hadron level before detector 
simulation, initial state radiation is turned off. In the case of the three-jet rate the detector 
corrections are Dh = 0.909 and Di = 0.843, with a statistical uncertainty at the few per mille 
level. 

The flavour tag can introduce a bias on the measured observables. The bias is estimated 
by computing the observable before and after applying the tag to the MC sample which passed 
the event selection cuts, for every flavour. As an example, Tb = 0.894 is found for the three-jet 
rate. Taking the product of ratios of detector corrections, tagging biases and purities the overall 
detector corrections amount to Df,/iTi,Vb = 0.801, Dc/iTcVc = 0.112 and TiVi = 0.045, again 
within a statistical accuracy of a few per mille. 



4.4 Hadronization corrections 

The perturbative predictions are corrected for hadronization effects by computing the relevant 
observables at parton and at hadron level, including final state photon radiation off quarks. 
Several Monte Carlo models based on the parton shower approach plus subsequent string or cluster 
fragmentation are employed for this purpose. For the nominal analysis, the same generator as for 



the full simulation is used, which is HVFL |22|. As already mentioned in Section [4.1| , this generator 
is based on the JETSET 7 . 4 parton shower model plus string fragmentation, and in addition the 
decay modes and lifetimes of heavy hadrons are adjusted to agree with recent measurements. 
Heavy quarks are fragmented using the Peterson et al. model. The fragmentation parameters for 



this generator are determined from a global fit to hadronic Z decay data as described in [|12[, using 
the present knowledge of heavy quark physics. The modelling of heavy quark physics represents 
an important part of the systematic uncertainty of the result. 

In the JETSET parton shower model mass effects are introduced by kinematic constraints to the 
phase space at each parton branching in the shower evolution as well as by the matching to the 
leading order matrix element at the first parton branching. The value of Mi, (pole mass) is set to 
4.75 GeV/ c^. Thirty million events without detector simulation were generated and analyzed. Such 
a large sample is necessary in order to reduce the statistical error of the hadronization corrections 
below the one per mille level. These corrections were computed for each flavour individually. 

The ratios of hadronization corrections, Hb/i, are listed in Table |^. The corrections are rather 
sizeable for almost all the observables; in most cases they are of the same size as or larger than 
the expected mass effect. Only the three-jet rate and the first two moments of the distribution 
have corrections at the 2% level or below. For all the other event-shape variables the corrections 
are of the order of 10% or even larger. It has been found that the deviation from unity is almost 
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Table 2: Results for -R^"* with statistical errors, including the uncertainty from the MC statistics. The results 
are obtained with the hadronization corrections as predicted by HVFL. Also given are the corrections H^°{'', 
computed from hadrons directly originating from the string, before any decays (7iod=no decays). 






ppcrt 
bd 




Tjnod 
^b/l 


R3 


0.974 ±0.005 


0.980 


0.987 




0.910 ±0.002 


1.134 


1.003 


Ci 


0.894 ±0.002 


1.169 


1.009 




0.955 ±0.005 


1.023 


0.987 




0.979 ±0.010 


0.984 


0.984 




0.831 ±0.001 


1.306 


1.025 


Bw2 


0.929 ±0.003 


1.088 


0.986 



entirely due to B hadron decays, which change the distributions mainly in the two-jet region. As 
can be observed from Table the same corrections, computed taking only hadrons which stem 
directly from the string before any decay, come close to unity within one or two percent. This is 
in agreement with expectations from recent calculations of nonperturbative power-law corrections 
to moments of event-shape distributions for massive quarks p3[. The strong dependence on the 



decays of heavy quarks is taken into account when estimating systematic uncertainties in the 
following. 



4.5 Measurement of R^^ 

The correction procedure defined in Eqn. ^ has been carried out individually for every data taking 
year. The hadronization corrections are taken from the HVFL generator and do not depend on 
the year. Good agreement between the values of the corrected observable B^i"^ in the different 
years is found, and therefore the final value is obtained from the average over years, weighted 
according to the statistical errors. The confidence level for the combination is 50%. In order 
to obtain the final values for B^'^^ , the contributions from anomalous triangle diagrams are first 
subtracted from the measured ratio B^i^^ , and then the error due to the finite MC statistics for the 
evaluation of purities and detector and tagging corrections is added to the statistical error of the 
data. This statistical error from MC is obtained by evaluating the scatter of results when repeating 
the analysis with a large number of MC subsamples of smaller size, and then extrapolating the 
standard deviation thus found to the actual size of the original MC sample. The results are listed 
in Table |. 

The dependence on the resolution parameter y^ut of the perturbative ratio B^'^^ for the three- 
jet rate is indicated in Fig. 0. The data are compared to the predictions of the parton shower 
models of PYTHIA 6.1, which is based on JETSET, and HERWIG 6. 1 (H. The MC models are in 
reasonable agreement with the measurement at large ?/cut values, where the statistical uncertainty 
is large, however. At intermediate resolution parameters these models predict a lower ratio than 
observed. In the lowest Ucnt region, where the resolution parameter approaches the scale of the h- 
mass, a turn-over is observed, which might be interpreted as a short-coming of the approximation 
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^ cut 

Figure 2: Comparison of the ?/cut dependence of the measured ratio i?}^"* for the three-jet rate to the predictions 
of parton shower models as well as next-to-leading order (NLO) perturbative QCD for two different values of the 
5-quark mass in the MS scheme. The measurement is obtained using HVFL for the hadronization corrections. The 
errors are statistical only. 

for mass effects as implemented in these models. 

In addition, in Fig. ^ the next-to-leading order perturbative QCD predictions for two different 
values of the running 6-quark mass in the MS scheme are given. The data clearly favour a mass 
close to 3 GeV/c^. 

4.6 Systematic uncertainties 
4.6.1 Detector and physics modelling 

Systematic uncertainties can arise from imperfections of the implementation of the physics 
processes as well as the description of the detector performance. An important physics parameter 
for the correct description of three-jet rates as well as tagging purities is the gluon splitting rate 
into heavy quark pairs (66, cc). The MC was reweighted in order to reproduce the most recent 
measured values and the experimental uncertainties were propagated to a systematic error on 

ppert 



ALEPH 




• Data, Stat errors only 
PYTHIA 6.1 



HERWIG6.1 

NLO, mi,(M^)=3 GeV/c^ 

NLO, mi,(M^)=5 GeV/c^ 
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The MC predictions for the background efficiencies depend on the mean number of VDET 
coordinates per track and on the assumed impact parameter resolution. It is observed that in 
the MC simulation more tracks are accepted for use by the 6-tag algorithm than in the data, and 
that the resolution function in impact parameter significance is not perfectly modelled. In order 
to improve the agreement between data and MC, tracks in the MC are deleted randomly, and 
a smearing algorithm is applied to the impact parameter significance of some of the MC tracks. 
A detailed description of this procedure is given in Ref. The systematic uncertainty related 
to this procedure is estimated by repeating the whole analysis once without the track deletion, 
and once without any smearing. Half of the observed deviations in the final result are taken as 
uncorrelated errors and added in quadrature. 

In the case of the three-jet rate additional cuts on the jet quality were studied. Removing the 
cut on the minimum jet energy results in a relative shift of —0.3% for . As a further check of 
the quality of the MC description of jet quantities, the requirement of at least three charged tracks 
per jet is applied. This causes a relative change in the result of —0.3%. A relative systematic 
uncertainty of 0.3% related to the modelling of three-jet events is therefore assumed. 

The relative uncertainties due to the physics and detector modelling as described above are 
listed in Table |^. Adding all these uncertainties in quadrature results in a final experimental 
systematic error (exp) as quoted in Table 11 for all the observables under consideration. 



4.6.2 Hadronization 

The uncertainties from the modelling of the hadronization are typically evaluated by computing 
the hadronization corrections with different MC generators. However, if these uncertainties are 
to be meaningful, it must be verified that the various models give a good overall description of 
hadronic Z decay data, and in particular of quantities relevant to the analysis. 

As has been shown in Section |4.4| , the B hadron decays have a large impact on the size of the 
hadronization corrections. Because in the tuning of the MC hadronic final states are analyzed 
after all decays, differences in the description of hadron decays can lead to differences in the 
tuned fragmentation parameters. In fact, comparing the predictions of the standard PYTHIA 6. 1 
and the HVFL generator, both based on the JETSET parton shower and string fragmentation and 
tuned to ALEPH data, differences in the fragmentation parameters are found, which translate to 
a variation in the hadronization corrections, even before considering decays. In order to assess an 
uncertainty related to the string fragmentation parameters and subsequent decays, the variation 
in the final result when using HVFL or PYTHIA for the hadronization corrections is taken as a 
systematic uncertainty. 

Another quantity relevant for this analysis is the h fragmentation function, which describes 
the energy fraction transferred to the B hadrons during the fragmentation process. In particular, 
the lower energy tail of the fragmentation function is important for three-jet events. Recently, 
new measurements of the h fragmentation function have been performed [p5| , |26| . Both HVFL 



and PYTHIA use the Peterson model in their standard setup. In order to test the sensitivity to 
the h fragmentation function, the measured distribution was included in a global tuning of 
the PYTHIA fragmentation parameters, and a new set of hadronization corrections was computed 
with those parameters. In Ref. |^ it is shown that the Lund fragmentation scheme combined 
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with a model by Bowler [17, 27] gives a good description of the b fragmentation function. The 
PYTHIA MC was therefore tuned to the data again, but now using the Bowler fragmentation scheme 
instead of the Peterson model, and new hadronization corrections were derived accordingly. The 
results obtained with the retuned parameter set and the Bowler fragmentation function were then 
compared to the standard PYTHIA predictions. The maximum deviation is taken as an estimator 
of the uncertainty related to the b fragmentation function. 

The HERWIG 6 . 1 MC generator was used for the study of a fragmentation model different 
from the string approach. As discussed above, differences in the description of hadron decays 
can effectively lead to differences in the pure fragmentation parameters, after tuning. In order to 
reduce the sensitivity to this effect, and to study purely the difference between string and cluster 
fragmentation, the variation in the hadronization corrections for light quarks (uds) only has been 
propagated into an uncertainty on the final result. 

The three components of uncertainty from the hadronization process as described above are 
listed in Table ^ They are then added in quadrature and quoted as the hadronization uncertainty 
(had) in Table ^. It is the dominant uncertainty for all variables. 



4.6.3 Additional systematic checks 

A number of additional systematic checks have been performed in order to evaluate the stability 
of the results. They are listed in Table |^. 

In Eqn. ^ the world average values for i?f, and Rc enter. Changing them within their errors 
has a negligible impact on the results. 

An additional factor which contributes to Eqn. ^ is -Rcf"^*, which is the ratio of the observable 
at parton level for c and light quark events. Because of the smaller c-quark mass this ratio is set 
to unity in the default analysis. Considering a c-quark mass of rric = 1.4GeV/c^ leads to very 
small variations. Only in the case of the total jet broadening is the observed shift larger than the 
statistical uncertainty; in that case the shift is therefore added in quadrature to the experimental 
uncertainty. 

A possible bias from the 6-quark mass in the parton level simulation is estimated by switching 
off the correction to the massive matrix element for the first branching of a massive quark in the 
parton shower. For most of the variables, and in particular for the three-jet rate, this results 
in small variations. As a further cross check, a variant of the JETSET fragmentation model was 
studied in which the parton shower cutoff parameter Qq is increased to 4 GeV. This setting leads 
to a parton multiplicity of 4 on average, which is close to the multiplicity used in the analytical 
next-to-leading order calculations. However, it represents a drastic change in the model, and such 
large cutoff values do not agree with the basic ideas of the parton shower approach. Nevertheless, 
recomputing Rf^^^ with the hadronization corrections obtained from this model leads to variations 
that are well within the fragmentation uncertainty quoted above. 

In order to check the stability of the result with respect to the chosen 6-tag working point, 
the cut on the output of the 6-tag algorithm was varied over a wide range, which resulted in 
changes of the purity Vb from 64% to 95%. For most of the variables, in particular the three-jet 
rate, the observed deviations in the final result are within the tracking uncertainty as described 
above. However, in those cases where the deviation exceeds this uncertainty, half of the maximum 
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Table 3: Variations (in percent) of R^^'^^ with respect to the nominal value for all systematic studies performed 
(PY=PYTHIA). The first block indicates the components of the statistical uncertainty from data and MC, in the 
second block the experimental uncertainties due to the physics and detector modelling are given, the third block 
lists the contributions to the systematic uncertainty due to the modelling of hadronization, and the last block gives 
the results of further systematic checks. 







Ti 














statistical uncertainties 


Stat, error Data 


0.44 


0.19 


0.15 


0.43 


0.90 


0.11 


0.28 


Stat, error MC 


0.28 


0.14 


0.11 


0.28 


0.55 


0.07 


0.20 




experimental uncertainties 


No hit smearing 


-0.47 


-0.26 


-0.22 


-0.42 


-0.50 


-0.14 


-0.36 


No track deletion 


0.71 


0.38 


0.29 


0.78 


1.10 


0.17 


0.57 


^ch/Jet > 


-0.29 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


iVch/Jet > 3 


-0.28 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


Gluon splitting 


-0.14 


0.10 


0.16 


-0.18 


-0.44 


0.25 


-0.03 




hadronization uncertainties 


n CTTl Tlfl TTI f^'tf^TQ. 
L idtilii. Ucxi CXliiC uCi o 
















and B decays 


-1.00 


0.93 


1.42 


-0.58 


-0.65 


3.20 


-0.16 


Fragm. model 


0.64 


-3.82 


-4.33 


0.67 


1.56 


4.42 


-0.89 


b fragm. function 


0.27 


-0.33 


-0.41 


0.14 


0.42 


-0.56 


-0.22 




additional systematic checks 


c-quark mass effect 


-0.03 


-0.12 


-0.15 


-0.09 


-0.02 


-0.34 


-0.11 


PY 6.1, 0{as) massless 


0.40 


-2.03 


-2.32 


0.14 


0.60 


-2.90 


-0.21 


PY 6.1, Qo = 4 GeV 


-0.86 


1.10 


1.47 


-0.82 


-0.66 


3.69 


0.27 


-log^o Puds > 1-2 


-0.25 


-0.57 


-0.71 


-0.05 


0.13 


-0.82 


-0.45 


- logio Puds > 4.0 


0.27 


0.82 


0.77 


0.69 


0.67 


0.67 


0.56 


Dmeas 
^bl 


0.03 


-0.04 


-0.11 


0.14 


0.30 


-0.14 


-0.09 



deviation under variation of the 6-tag working point is taken as a systematic uncertainty due to 
the tracking. 

An attempt was made to verify with data the purity Vb determined from MC. The approach 



is based on the double tag method as described in It is possible to extract the hemisphere 
tagging efficiencies, which then can be related to a global event tagging efficiency if the hemisphere 
correlations are known. The same working point as in the nominal analysis was chosen. In 
a first check, Rf, and the hemisphere 6-tag efficiency are measured from the data, whereas 
all the background efficiencies and the hemisphere correlations are taken from MC. The results 
are Rb = 0.2190 ± 0.0013(stat), which is consistent with the world average value, and = 
0.5380 ±0.0024(stat). Next Rb is fixed to the world average value, and as well as the hemisphere 
correlation pb are measured. The result can be translated to the global efficiency by the relation 
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Table 4: Results for R^^^^ with statistical (stat), experimental (exp) and hadronization (had) uncertainties. 



o 


Tjpcrt 
R-bd 


±(stat) 


±(exp) 


±(nad) 


-n-S 


n Q7zL 




yj.yjyjo 




Ti 


0.910 


0.002 


0.004 


0.036 


Ci 


0.894 


0.002 


0.004 


0.041 


l/3i 


0.955 


0.005 


0.005 


0.009 


1/32 


0.979 


0.010 


0.007 


0.017 




0.831 


0.001 


0.005 


0.046 




0.929 


0.003 


0.003 


0.009 



efc = 2e^ — {e^y (1 + pb), and then the purity is obtained as 



Rb€b + Rc^c + (1 — -Rfc — Rc) ^uds 

The result is Vb = 0.8310 ±0.0007(stat), which is in very good agreement with the value extracted 
from the MC, V^^ = 0.8305. 

In order to test the overall correction method, an independent set of MC events was analyzed 
in the same way as the data. Within statistical errors, the extracted -R^f'^* values are in agreement 
with the "true" values as found at the parton level of the MC prediction. 

Finally, the whole measurement is repeated using the method of Ref. ||^, where R^'^'^^ is 
measured instead of Rf^j^^- Both methods have the same statistical accuracy, but the latter is 
more simple from the point of view of systematic uncertainties. For the purpose of measuring 
^meas ^ must also be applied. This is achieved by the requirement — log^Q Puds 

< 0.5, 

which selects light quark events with an efficiency of 66% and a purity of Vi = 88%. The correction 
procedure and the correction factors needed are different; nevertheless the results obtained differ 
from the nominal ones by at most 0.3%. 



4.7 Results for R^J* 



As stated in Section [4.4| , the nominal results are obtained when using HVFL for the hadronization 
corrections. The results for all observables can be found in Table ^, together with statistical 
and systematic uncertainties. The result for the three-jet rate is in agreement with that of 
Ref. 0. Most of the event-shape variables have a better statistical accuracy than the three-jet 
rate. However, the hadronization uncertainties are rather large. 

For the extraction of the 6-quark mass it is required that the uncertainty on the measured ratio 
be at the 1% level, and that at the same time the size of the hadronization correction (Table be 
smaller than the actually measured mass effect, given by the deviation of R^^^^ from unity. These 
requirements leave only the first moment of the ys distribution and the three-jet rate for further 
investigation. 
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5 Determination of the 6-quark mass 



5.1 Theoretical predictions 



The NLO prediction for -R^^'^* as a function of the observable O and the 6-quark mass can be cast 
in the following form: 



i?r(0) = l+r,(/i) 



6o(r,(/i),0) + ^6i(r,(/i),0) 



1 + K(rb(^)) 



(6) 



where rb(yu) = m^(/i)/M|, and mf,{fi) is the running 6-quark mass as defined in the MS 
renormalization scheme at the renormalization scale fi. The ratio can also be expressed in terms 
of the pole mass Mf,, 



b ) 



0) + ^6f(rf,0) 

/TT 



(7) 



with 



M^/M^. The two predictions are equivalent at this order. The coefficient functions 



6o,i for the two schemes can be related to each other by expressing the pole mass in terms of the 
running mass, i.e., 



^ + ^ (t-^ ^'"^ 



The coefficient functions for the three-jet rate were computed in Ref. |I|]. For all the other variables 
they were obtained using the MC generators ZBB4 



28] and EVENT pg. ZBB4 allows for the 
integration of the fully differential NLO matrix elements including mass effects, whereas EVENT 
contains the massless expressions The latter has been extensively employed at LEP for the 
calculations of the NLO predictions of event shapes used for the measurements of the strong 
coupling constant. The differential cross sections for any infrared and coUinear safe observable O 
for either b or d primary quarks in e"'"e~ annihilation can be written as 



1 da, 



b,d 



<4 dO 



Zn VzTT/ 



(9) 



Here o"^°j is the total cross section for either b or d quark production. The functions Ab, Bh are 
computed with ZBB4 for fixed values of the 6-quark mass in the pole mass scheme, whereas A^ 
and Bd are obtained from EVENT. Taking the ratio of the cross sections for b and d quarks and 
expanding up to NLO, the following relationships can be found: 



^b % 



Ad 



- 1 



'b 



BbAd — BdAb 



(10) 



:ii^ 



In order to get the functional dependence of b^i on the 6-quark mass, the coefficients Ab and 
Bb were first computed in very high statistics runs for values of Mb = 1, 2, 3, 4, 5, 6 GeV/c^. 
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Figure 3: Parametrisations of the 6-quark mass dependence for the LO (top) and NLO (bottom) contributions 
to R^^^^ in the pole mass (left) and running mass schemes (right) for the first moment of the ys distribution. The 
points indicate the result of the MC integration, and the full line the parametrisations. The NLO contributions 
are evaluated using as{M-z) = 0.119. 



The expressions (^U|) and (|TT]) were then fitted over the relevant range of 3 to 6 GeV/c^ using 
parametrisations of the type ci + C2r^ + C3 Inr^. This functional form follows closely that used in 
. It has been checked that different parametrisations lead to very small changes in the extracted 
mass values. 

In Fig. ^ an example of these fits is shown for the first moment of the distribution. With 
these parametrisations it is possible to estimate the actual quark mass effects in leading order and 
its next-to-leading order corrections. In Table |^ a list of these LO and NLO contributions is given 
for all variables for both the running and the pole mass schemes. 

Possible evaluations of uncertainties in the mass extraction because of the limited accuracy 
of the NLO predictions are discussed next. Missing higher order corrections are estimated by 
extracting firstly the pole mass from the perturbative expression in the pole mass scheme, then 
translating that result into a running quark mass mb{'mb) at the 6-mass scale, and finally running 
this mass up to the Mz scale, using the relation 



where = a.(Mz)/[l + a,(Mz)/3o ln(/iVM|)/(47r)], /3o = 11 - 2n//3, % = 5 and 70 = 2. 

The world average value ^^(Mz) = 0.119 ± 0.003 |^ is used for the strong coupling constant. 
At NLO the two methods of mass extraction are equivalent. However, the sensitivity to higher 
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Figure 4: Renormalization scale dependence for the extracted 6-quark mass in the running and pole mass scheme, 
for the observables i?3 (left) and y^^ (right). The error bar spans the range of the scheme uncertainty. The full line 
indicates the renormalization scale dependence of the average computed from the running and pole mass schemes. 



order contributions is different for tlie two scliemes, wliicli can tlien lead to different results. The 
shift in the mass when using the pole mass scheme is — 0.14GeV/c^ and — 0.4GeV/c^ for the 
distribution and the three-jet rate, respectively. 

In addition, the effects of uncalculated higher order terms can be estimated by a change in 
the renormalization scale. For the central value /i = Mz is employed, and the systematic error is 
taken to be half of the range of mass values (average of running and pole mass schemes) found 
when varying /i from O.lMz to 2Mz. This results in a scale error of 0.15 GeV/c^ and 0.06 GeV/c^ 
for the j/3 distribution and the three-jet rate, respectively. The behaviour of the extracted mass 
in the two schemes as a function of the renormalization scale /x is illustrated in Fig. ^. 

The uncertainty on the world average value for the strong coupling constant has a negligible 
impact, about 0.01 GeV/c^, on the measured 6-quark mass. 



5.2 Results for mb(Mz) 

Based on the predictions obtained above, the running 6-quark mass is determined from the 
measured ratio B^'^ for two observables, the first moment of the distribution and the three-jet 
rate. The results are listed in Table ^ together with the statistical and systematic uncertainties, 
which have been propagated from the corresponding uncertainties on I^'^^ ■ The statistical 
correlation between the two measurements is 0.66, while the hadronization uncertainties are almost 
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Table 5: Measured 6-quark mass mb(Afz) in the MS scheme with statistical (stat), experimental (exp) and 
hadronization (had) uncertainties. In the last column the measured pole mass Mjj is listed. 






mfe(Mz) [GeV/c2] 


±(stat) 


±(exp) 


±(had) 


Mb [GeV/c2] 




2.76 


0.28 


0.28 


0.62 


3.65 




3.34 


0.22 


0.22 


0.38 


4.73 



100% correlated. A similar difference of measured masses was reported in Ref. [Q, where different 
clustering algorithms for the three-jet rate were studied. The results obtained for the pole mass 
by fitting the predictions in the pole mass scheme are also given in Table ^. 

As is observed from Table the first moment of the distribution allows for a 
mass measurement with better statistical and experimental accuracy than the three-jet rate. 
Furthermore, the hadronization uncertainty is considerably smaller. Therefore the result from 
this observable is quoted as the final 6-quark mass value. The theoretical uncertainty is estimated 
as described in Section |5.1| by evaluating the impact of the uncertainty on the strong coupling 
constant and the renormalization scale variation. Because the pole and the running mass scheme 
are equivalent at NLO, the average of the values found for the two schemes is quoted as the final 
result, and half of the difference is taken as an additional theoretical systematic uncertainty due 
to the scheme ambiguity. This leads to a measurement of the 6-quark mass of 

mb{Mz) = [3.27 ± 0.22(stat) ± 0.22(exp) ± 0.38(had) ± 0.16(theo)] GeV/c^ . 

The result obtained is in agreement with the measurements of DELPHI |^ and Brandenburg et 
al. p. The world average value from low-energy measurements is mb{mb) = 4.2 ± O.lGeV/c^ 



3^. Using a two- loop running equation for the quark mass, this average value can be translated 
to mfe(Mz) = 2.89 ± 0.08GeV/c^, which is also in agreement with the result obtained here. This 
comparison is illustrated in Fig. ^ together with the other results. The result for the pole mass is 

Mb = [4.73 ± 0.29(stat) ± 0.29(exp) ± 0.49(had) ± 0.18(theo)] GeV/c^. 

In this case the theoretical uncertainty is estimated from the renormalization scale variation only. 



6 Test of the flavour independence of cXg 

All the previous considerations are based on the assumption of flavour independence of the strong 
coupling constant, which in particular means that the strong coupling is the same for b quarks as 
for light quarks (/ = uds), i.e., = a[. However, the measurement of -R^^"^* can also be used as a 
test of the flavour independence, if the 6-quark mass is known. The ratio of the strong coupling 
constant for b and light quark events can be written up to NLO as 

^ = [RlT - Kinm))] + a,^^ [Rlf - KiniMz)) - 1] , (13) 
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Figure 5: Comparison of the ALEPH result for mb{Mz) with the world average value of low-energy measurements 
for mb{mb), which is evolved up to the A/z scale using a two-loop evolution equation with as(Mz) = 0.119 ± 0.003. 
Also shown are the measurements by DELPHI |^ and Brandenburg et al. The inner error bars indicate the 
quadratic sum of the statistical and experimental uncertainties. The three points at the Z pole are separated for 
clarity. 



where K{rb{Mz)) is the mass correction in Eqn. |^ and ai = 6.073 when using the first moment 
of the 1/3 distribution as observable. Equation |1^ is derived assuming = a[{l + 6), 6 <^ 1, 
and neglecting all terms of order (a^)^ and a[6, so a<j(Mz) = 0.119 ± 0.003 can be used for 
the coefficient of the next-to-leading order term. The coefficient ai was computed with EVENT. 
Taking an average 6-quark mass of mfe(mb) = 4.2 ± O.lGeV/c^ and taking into account 



the scale uncertainties and scheme ambiguities as described in Section pTT] , the term K{rb{Mz)) 



amounts to —0.041 ± 0.003. Inserting the measured value R^^^^ for the observable y^-^ results in 

^ = 0.997 ± 0.004 (stat) ± 0.004 (exp) ± 0.007 (had) ± 0.003 (theo) , 
a. 

which is a confirmation of the flavour independence at a precision level of 1%. This constitutes 
an improvement in precision by a factor of two compared to a previous ALEPH analysis ||33| . 
Recently, similar results were obtained by other experiments H, 0] • 
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7 Conclusions 



The effect of the 6-quark mass has been studied on the ratios of observables in b and hght quark 
decays of the Z boson in e"'"e~ annihilation. Taking the first moment of the distribution, which 
has the smallest hadronization corrections and systematic uncertainties, a 6-quark mass of 

mfc(Mz) = [3.27 ± 0.22(stat) ± 0.22(exp) ± 0.38(had) ± 0.16(theo)] GeV/c^ 

is found, which is in agreement with recent measurements by other experiments [|, H and 
with measurements at lower energies. The three-jet rate defined by the DURHAM algorithm with 
Vcut = 0.02 turns out to be more sensitive to the details of the hadronization modelling. 

The determination of the ratio of the first moment of the distribution in b and light quark 
events has alternatively been employed, along with previous measurements of the 6-quark mass, 
to test the flavour independence of the strong coupling constant with a precision of 1.0%. 
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